***************************************************
****----			Table 1	  			   ----****
***************************************************

clear all
cap log close 
set more off

****---- Change directory

global data "/Users/Replication material/Data"

****---- Analysis

use "$data/dataset_replication.dta", clear

global controls Capoluogodiprovincia  population_density_1911 expenditures_1912_PC surplus_tot_1912_PC exp_police_hygiene_1912_PC exp_just_secur_1912_PC  exp_edu_1912_PC exp_public_work_1912_PC share_literate share_literate_male share_literate_female share_male


qui reghdfe rate_01_11 aug_dec_18_sd, cluster(distretto_FE) absorb( i.Zonaaltimetrica  i.N_contribuenti_quartile##i.distretto_FE)
est store reg11
qui reghdfe population_density_1911 aug_dec_18_sd, cluster(distretto_FE) absorb( i.Zonaaltimetrica  i.N_contribuenti_quartile##i.distretto_FE)
est store reg12
qui reghdfe Zonaaltimetrica aug_dec_18_sd, cluster(distretto_FE) absorb(  i.N_contribuenti_quartile##i.distretto_FE)
est store reg13
qui reghdfe N_contribuenti_quartile aug_dec_18_sd, cluster(distretto_FE) absorb( i.Zonaaltimetrica i.distretto_FE)
est store reg14
qui reghdfe Capoluogodiprovincia aug_dec_18_sd, cluster(distretto_FE) absorb( i.Zonaaltimetrica  i.N_contribuenti_quartile##i.distretto_FE)
est store reg15

qui reghdfe expenditures_1912_PC aug_dec_18_sd, cluster(distretto_FE) absorb( i.Zonaaltimetrica  i.N_contribuenti_quartile##i.distretto_FE)
est store reg21
qui reghdfe surplus_tot_1912_PC aug_dec_18_sd, cluster(distretto_FE) absorb( i.Zonaaltimetrica  i.N_contribuenti_quartile##i.distretto_FE)
est store reg22
qui reghdfe exp_police_hygiene_1912_PC aug_dec_18_sd, cluster(distretto_FE) absorb( i.Zonaaltimetrica  i.N_contribuenti_quartile##i.distretto_FE)
est store reg23
qui reghdfe exp_just_secur_1912_PC aug_dec_18_sd, cluster(distretto_FE) absorb( i.Zonaaltimetrica  i.N_contribuenti_quartile##i.distretto_FE)
est store reg24
qui reghdfe exp_edu_1912_PC aug_dec_18_sd, cluster(distretto_FE) absorb( i.Zonaaltimetrica  i.N_contribuenti_quartile##i.distretto_FE)
est store reg25

qui reghdfe exp_public_work_1912_PC aug_dec_18_sd, cluster(distretto_FE) absorb( i.Zonaaltimetrica  i.N_contribuenti_quartile##i.distretto_FE)
est store reg31
qui reghdfe share_literate aug_dec_18_sd, cluster(distretto_FE) absorb( i.Zonaaltimetrica  i.N_contribuenti_quartile##i.distretto_FE)
est store reg32
qui reghdfe share_literate_male aug_dec_18_sd, cluster(distretto_FE) absorb( i.Zonaaltimetrica  i.N_contribuenti_quartile##i.distretto_FE)
est store reg33
qui reghdfe share_literate_female aug_dec_18_sd, cluster(distretto_FE) absorb( i.Zonaaltimetrica  i.N_contribuenti_quartile##i.distretto_FE)
est store reg34
qui reghdfe share_male aug_dec_18_sd, cluster(distretto_FE) absorb( i.Zonaaltimetrica  i.N_contribuenti_quartile##i.distretto_FE)
est store reg35
           
esttab reg11 reg12 reg13 reg14 reg15  , keep(aug_dec_18_sd) se r2 staraux collabels(, none) ml(,none) cells(b(star fmt (%9.3f)) se(par)) star(* 0.10 ** 0.05 *** 0.01) replace
esttab reg21 reg22 reg23 reg24 reg25  , keep(aug_dec_18_sd) se r2 staraux collabels(, none) ml(,none) cells(b(star fmt (%9.3f)) se(par)) star(* 0.10 ** 0.05 *** 0.01) replace
esttab reg31 reg32 reg33 reg34 reg35, keep(aug_dec_18_sd) se r2 staraux collabels(, none) ml(,none) cells(b(star fmt (%9.3f)) se(par)) star(* 0.10 ** 0.05 *** 0.01) replace



